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O ' Abstract 

o 

^2 ■ The lectures presented by one of us (IP) at the Les Houches summer school dealt 

with the scaling properties of high Reynolds number turbulence in fluid flows. 

G\ ■ 

"^5 ■ The results presented are available in the literature and there is no real need to 

£>y reproduce them here. Quite on the contrary, some of the basic tools of the field 

and theoretical techniques are not available in a pedagogical format, and it seems 
§ ' worthwhile to present them here for the benefit of the interested student. We 

begin with a detailed exposition of the naive perturbation theory for the ensemble 
averages of hydrodynamic observables (the mean velocity, the response functions 
and the correlation functions) . The effective expansion parameter in such a theory 
is the Reynolds number (Re) ; one needs therefore to perform infinite resummations 
to change the effective expansion parameter. We present in detail the Dyson- Wyld 
line resummation which allows one to dress the propagators, and to change the 
effective expansion parameter from Re to 0(1). Next we develop the "dressed 
vertex" representation of the diagrammatic series. Lastly we discuss in full detail 
the path-integral formulation of the statistical theory of turbulence, and show 
that it is equivalent order by order to the Dyson- Wyld theory. On the basis of the 
material presented here one can proceed smoothly to read the recent developments 
in this field. 
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1. INTRODUCTION 



The exact description of fluid flows is provided by the specification of the velocity field 
u(r, t) as a function of space and time. Such a description is natural for "laminar" flows of 
fluids whose velocity field is not strongly fluctuating in space and time. On the other hand, 
for "turbulent" flows, which are states of fluids whose velocity field u(r, t) appears highly 
erratic as a function of r and t, it is more sensible to consider a statistical description. 
The statistical description is particularly useful in situations where the flow is driven by 
stationary external agents. Examples of such stationary driving mechanisms are found in 
shear turbulence, grid turbulence, convective turbulence, etc. It is believed that in such 
flows it makes sense to consider long time averages of functions of the velocity field u(r, t). 
For example, if the driving is arranged such that the fluid as a whole is confined in some 
volume V, the mean velocity must vanish: 



X V ' " V J T^oo2TJ-T v ' ' v ' ' 

The pointed brackets indicate here an average with respect to space and time. Likewise, 
correlation functions can be defined: 



If the state of the fluid is statistically homogenous in space and time, the averaging process 



One of the most important open problems of statistical physics is whether the statistical 
description of turbulence displays universal characteristics. Many important concepts were 
developed in the attempt to answer this question. It was recognized by Richardson in 
the '20's that quantities which depend on one space point have no chance to claim univer- 
sality. The velocity field u(r, t) and all its functions are fatally sensitive to the large scale 
components which are introduced by the external driving agents. In this picture the driving 
acts on the large scales, leading to maximally large velocity fluctuations on the macroscale L, 
which is known as the "outer" scale of turbulence. The typical scale of the velocity fluctua- 
tions is denoted as Ul- The order of magnitude of u(r, t) at every r is about Ul- Richardson 
suggested that the quantities that may have universal aspects are the velocity differences 
8u(r + R, r, t) = u(r + R, t) — u(r, t). In taking differences one eliminates the contribution 
of the large scale motion, and may hope to see universal properties that are independent of 




(1.1) 



(u(r + r',t + r)u(r,t)) 



(1.2) 




( |1.2|) results in a function of r' and r only. 
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the driving mechanism. It was one of the main suggestions of Kolmogorov that for R that 
is much smaller than L and much larger than a "dissipative" scale rj, statistical averages of 
functions of 8u(r + R, r, t) may show universality. The dissipative scale rj depends on the 
kinematic viscosity of the fluid z/, and on the single most important dimensionless number 
in the theory of turbulence, i.e. the Reynolds number Re. This number is defined as || 

i?e = ^, (1.3) 

and in the Kolmogorov approach it is asserted that r//L ~ i?e~ 3 / 4 . In other words, as 
Re— > oo there is a larger and larger range of scales R for which functions of 5u(r + R, r) are 
expected to show universal properties. 

In most experimental applications it is customary to consider the longitudinal component 
of 5u(r + R, r, t) only: 

R 

5u(r + R, r, t) = 5u(r + R, r, t) ■ — . (1.4) 

R 

The essence of the Kolmogorov approach to the issue of universality in turbulence, which 
was suggested in 1941 Q (K41) is that it is possible to construct a local theory with one 
universal scaling exponent. This universal exponent was ascribed to 5u(r + R, r), stating 
that 

5u(r + R, r) ~ R 1/3 . (1.5) 

It was asserted that 5u (r + R, r) is the only scaling field that has to be considered, in the 
sense that the structure functions S n (R) which can be formed from 5u(r + R, r, t) satisfy 
the scaling laws 

S n {R) = (([6u{T + R,r,t)] n )) 

~ (eR) nCn ~ (eR) n/3 (1.6) 

for values of R in the "inertial range" L ^> R ^> rj. In ( |1.6| ) e is the mean (really in the 
double bracketed sense of ( |1 . 1|) but defined by an overbar to save notation) of the dissipation 
field e(r, t), 

e(r, t) = ~ [d a u p (r, t) + d p u a (r, t)f . (1.7) 

This suggestion of Kolmogorov was immediately attacked both by theorists (e.g. Landau) 
and experimentalists (for a recent review see e.g. ||). Indeed, it is rather astonishing 
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that a problem like fluid turbulence, which suffers from very large fluctuations and strong 
correlations, should be amenable to such a simple description; even Kolmogorov himself 
revised his thinking and changed (|1.6|) to a more complicated form (which fell under attack 
as well). Indeed, one measurement that raised a lot of objections to the K41 approach is 
the measurement of the correlation function of the dissipation field K ee (r) 

K ££ (R) = (e(r + R, t)e(r, t)) ~ iT" . (1.8) 

where e(r, t) = e(r,t) — e. It was found in experiments that K ee (R) decays very slowly in 
the inertial range, with fi having a numerical value in the range 0.20-0.25 ||. It was claimed 
that the K41 theory required fi to vanish. According ly, there have been many attempts to 
construct models of turbulence to take (|1.8|) into account and to explain how the measured 
deviations (£„ — 1/3) in the exponents of the structure functions were related to fi. 

Our aim is to present a Navier-Stokes based theory of turbulence. This theory attempts 
to study the universality of of turbulence on scales belonging to the bulk of the inertial 
range using the equations of fluid mechanics. We are not going to review the variety of 
phenomenological ad hoc models which have proliferated in recent years. It is our belief 
that recent progress in the analytic theory brings us a long way towards understanding the 
scaling properties of turbulence, and it is important at this point to examine carefully what 
is known, to separate false starts from genuine progress, and to delineate the domain of firm 
results. In doing so we will also attempt to carefully assess the main open problems that 
await further research. 

The equations of fluid mechanics are partial differential equations, and it is natural to 
consider them as a classical field theory, and to attempt to apply techniques of perturbation 
theory and renormalization that found spectacular success in the context of quantum field 
theory. Indeed, such attempts were made with the help of various closure procedures like 
Kraichnan's Direct Interaction Approximation and many others (for review of different 
closure schemes see, e.g. |§). A more systematic approach was suggested to be found 
within diagrammatic perturbation approaches for non-equilibrium processes like the Wyld 
diagrammatic technique |J for the Navier-Stokes equation. Generalizations of the Wyld 
diagrammatic technique were suggested by Martin, Siggia, Rose [[TIJ and by Zakharov and 
L'vov 0. 

In performing such calculations one faces several serious difficulties. Firstly, turbulence 
has no small parameter. The naive perturbation expansion of the Navier-Stokes equations 
results in increasing powers of the Reynolds number, and this is rather inconvenient in a 
situation where one is interested in the limit Re3> 1. Wyld [f| found a way to overcome this 
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difficulty; by performing line renormalization like Dyson's quantum field analog one finds a 
perturbative theory in which the coupling parameter is of 0(1). This is an important step, 
but it does not save Wyld's theory from severe problems. The main problem of the Wyld 
theory is that it is written in terms of the correlation functions of the field u(r, t) itself. 
As mentioned above, such quantities are not universal, and they are dominated by large 
scale contributions. In addition, functions of u(r, t) itself are not Galilean invariant. Wyld 
attempted therefore to work in terms of Fourier transforms of these quantities. Indeed, 
the correlation functions in k, uj representation are expected to be universal in the regime 
1/f] ^> k ^> 1/L. Unfortunately, the theory in k, u representation involves integrals over the 
whole k, u range, and the theory picks up infra-red divergences that cannot be eliminated 
easily. 

Kraichnan has suggested that a theory that does not suffer from the problems of Wyld's 
expansion can be formulated in Lagrangian coordinates , F3 . From the point of view of 
developing the theory in terms of universal quantities this suggestion was correct. Unfortu- 
nately the perturbation theory which is based on the Lagrangian representation is not of the 
diagrammatic type, and higher order terms with its coefficients cannot be found simply on 
the basis of topological features of the diagrams. We shall not be interested here in theories 
that contain arbitrary truncations. Such theories are uncontrolled, and it is our hope that 
one can achieve a fully controlled theory of turbulence. We seek therefore a formulation that 
lends itself to a diagrammatic description and which comes as close as possible to using as 
building blocks functions of Su(r + R, r). A theory that has all the right ingredients was 
formulated by Belinicher and L'vov |I4| (see also [|HJ). This theory went under the name 
"quasi-Lagrangian" ; this name turned out to be unfortunate, since it left the impression of 



an approximate theory. Recently |16] we revisited this theory and reiterated that it is not 
an approximate theory, but rather an exact renormalization scheme that differs from the 
Wyld technique. 

The aim of these notes is to review the basic tools of the field theoretic diagrammatic 
approach to the Navier-Stokes problem. There are essentially no new developments in these 
notes, but they are meant to bring the student to a level of expertise that would allow a 
smooth entry into the latest developments. We develop in full detail the naive perturbation 
theory (sec. 2) for the Green's function and the 2-point correlation function. This theory 
is badly diverging in powers of Re. In sec. 3 we discuss the line resummation which results 
in the Dyson equation for the Green's function and the Wyld equation for the correlator. 
In these equations, which are also presented as infinite diagrammatic series expansions, the 
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effective coupling constant is of 0(1). We discuss briefly 3-point and higher order correlation 
functions, and in Sect. 4 we turn to the functional integral formulation of the renormalized 
perturbation theory. We show in detail that the diagrammatic expansion obtained in this 
formulation is identical, order by order, with the Wyld formulation, and that either one 
can be used at will JTU| , p!7|P^| . The interested student is invited then to continue reading 
Refs. fln||r9| , pD| in which the Belinicher-L'vov renormalization scheme is used to develop a 
consistent theory of normal and anomalous scaling in turbulence. 

2. NAIVE PERTURBATION THEORY 
A. Introduction 

As was mentioned in section [I], we are mostly interested in the calculation of averages 
and correlation functions, which are obtained in the long time limit of the forced Navier- 
Stokes equations. Since the equations contain a viscous damping, it is assumed that the 
balance between forcing and damping leads after a time to stationary statistics which are 
independent of the initial conditions. Roughly speaking the picture is of a nonlinear problem 
that generates its own ergodic measure, and the forcing is chosen such that the ergodic 
measure (at least as far as the calculated quantities are concerned) is a property of the 
Navier-Stokes equations and not of the forcing. In other words, the solution roams forever 
on a strange attractor, and the forcing is used mainly to keep the system away from the 
(stable) zero solution, but without ruining too much the measure on the attractor. 

In our formal scheme we will expand ergodic averages for the nonlinear problem in terms 
of averages of solutions of the linear problem (i.e with the nonlinearity discarded). Not 
surprisingly, since the measure of the nonlinear problem is very different from the measure of 
the linear problem, we will see that our expansion has very miserable convergence properties. 
Actually the naive expansion is diverges, and it cannot be truncated at any order. Its 
has meaning only for the series as a whole. Nevertheless, one can perform infinite partial 
resummations of this series to achieve a new series that has somewhat better properties, 
i.e. the terms remain all of the same order. We stress right at the beginning that many of 
our steps have dubious mathematical validity, and the fact that they are common practice 
in many field theories cannot be really taken as an excuse. At least we will not hide the 
difficulties. 

The Eulerian velocity field is denoted in our discussion as u(r, t). In later notation we 
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shall sometimes use the 4- vector notation u(x), where x = (r, t). In the mathematical 
literature it is customary to write the Navier-Stokes equations as an initial value problem 
without any type of forcing: 

7? + (u- V)u-zA7 2 u- Vp = 0, V-u = 0. (2.1) 

For the physicist the equations of fluid mechanics are a projection onto a set of slow variables 
of the many-body microscopic dynamics. In such a projection, when successful, the host 
of irrelevant degrees of freedom remain in the form of a "noise" term in the equation of 
motion of the slow variables. The basic assertion is that the irrelevant degrees of freedom 
are much faster, and reach local equilibrium on a short time scale which is irrelevant for the 
hydrodynamic description. Thus, the physicist writes Eq. (|2.1| ) in the form 

du/dt + (u- V)u- z/V 2 u- Vp = f, (2.2) 

where the fluctutation dissipation theorem dictates that the two-point correlation function 
of the "random" force f is related to the viscosity and the temperature of the fluid: 

f a (T,t)fp(r>,t>) = 2uT5(r - r')5(t - t')5 a/3 . (2.3) 

An overbar denotes an average with respect to the thermodynamic equilibrium ensemble 
with temperature T. 

In natural conditions turbulence arises either due to the instabilitiy of laminar flows (like 
in channel flow) or due to time dependent stirring (such as in a glass of water stirred by 
a spoon). For the sake of theoretical simplicity we will model all these situations with a 
stirring force 0(r, t). The properties of the turbulent flow may depend on the nature of the 
stirring force 4>(r,t). Since in experimental situations one finds that the statistics of the 
turbulent velocity field on the large scales are close to that of a Gaussian random ensemble, 
it is customary to take f(r,t) to be a Gaussian random force with spectral support in the 
small wave vectors and small frequencies. The properties of the correlation function of 
0(r, t) are best stated in k, u space: it is concentrated in the small k, uj region, i.e. k < 1/L, 
uj < U L where U L is the characteristic magnitude of turbulent velocity fluctuations. In order 
to model properly the generic properties of turbulent fluids one needs to assume that the 
correlation decays sufficently quickly to zero for k ^> 1/L, uj 3> U L / L. In r, t space it means 
that D(R,t), 

D a p(R,r) = ((<f> a (r + R,t + T)<l>p(T,t))) (2.4) 
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is constant for r ^ L, r ^ L/U and for higher values of R and r it decays quickly. 
One of the important questions that the theory should address is whether the properties of 
turbulence on scales much smaller than L are independent of the precise choice of D(R, r). 
Also, the possible effect of non-Gaussianity of the forcing should be understood. 
Adding the stirring force to the equation of motion we write 

du/dt + (u- V)u-z/V 2 u- Vp = f + . (2.5) 

The combination f + <f> will be referred to, when convenient, as f . 

In this work we deal with incompressible turbulence only, in which V • u = 0. We shall 
therefore project out from Eq. fl2.5| ) any longitudinal components. This is done with the 
help of the projection operator P which is formally written as P = — V~ 2 V x Vx . In tensor 
notation this operator is P a p = 5 a p — V _2 V a V / 3. Because of the inverse Laplace operator 
which is nonlocal in space, the application of P to any given vector field a(r) is non local: 

[Pa(r)] Q = J dv'P aP {Y - v')a p {r') , (2.6) 

where P Q/ g(r — r') is the inverse Fourier transform of P aj g(k) 

r r/k 

Papir - r') = J 7^p ex PH( r " r ') • k]^(k). (2.7) 
The tensor P Q/ g(k) in ^-representation is 

Pap(k) = 5 aP - ^k a k p (2.8) 

and in r-representation 

^-^R- i ' r, M."''' (2 ' 9) 

The operator P has the following properties: (i) V ■ P = P • V = 0. (ii) For any divergence- 
free vector b, Pb = b. Applying P to Eq. fl2.5p we find 

(d/dt-uV 2 )u + P(u- V)u = Pf . (2.10) 

Introduce the bare Green's operator Go via the relationship 



Using Go we can rewrite Eq. ( |2.10|) in the form: 
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u(r,i) = u (r,t) - zG (u • V)u , (2.12) 

where 

Uo(r,*) = -zG f . (2.13) 

The meaning of the application of Go to any vector field c(r, t) is understood by 

G c(r, t) = G%(r -r',t- t') * c(r', t) (2.14) 
= J dv'dt'Glpiv-T^t-t'ypiv'.t'), 

where the kernel Cr^g(r — r', t — t') is the Green's function of the linear part of Eq. ( p,10| ), 
and we have introduced the * operation according to the RHS of Eq. ( |2.14j ). The explicit 
representaion of — r',t — t') is given as the inverse transform in k and uj 

G° a/3 (v-v',t-t') (2.15) 
dkduo 



(2k) 

The function G^^k, uj) is 

The convention chosen is that G°^(k, u) is analytic in the upper half of the plane of complex 
uj. Correspondingly ^^(r — r',t — t') is zero for t' < t. This property is known as causality, 
and it stems from the phyiscal constraint that a response cannot come before its cause. 



It is obvious that the solution ( 2.12 ) involves integrals over all space in r, t representation 
because of the inverse operators in Go which is a factor function in k, uj space. The solution 
involves integrals also in k, uj space due to the nonlinear term which is a local operator in 
r, t space. Thus there is no formal advantage to either representation, and we shall write 
Eq. (|2.12|) in the representation invariant form 

u = u Q + ^G°*r-. (2.17) 

The * and the • operations depend on the representation. In r, t representation the 
* operation was defined by Eq. ( [2.14j ), and E^ means 2(u ■ V)u. For later applications we 
need to consider the operation EJJ in which different fields u and u are involved. We will 
interpret E£ in a symmetric fashion, i.e. 

r;g = -i[(u- V)u+(u- V)u] . (2.18) 



s 



In Fourier representation we write 



and 



where 



u(k, uj) = J drdt exp[— i(r ■ k + o;t)]u(r, t) , 



G° * c 



G i0(k,u)ci3(k,uj) 



rr u i - / rfk ' duj" 



r a /3 7 (g, q', q") = (27r) 4 (fc /3 5 a7 + fc 7 5 a/3 ) 

x <J(k + k' + k")5{u + uj' + uj' 1 



(2.19) 
(2.20) 



(2.21) 



(2.22) 



In these equations and below we introduced the four dimensional vector q = {k, uj}. The 
symmetry of ( [2.21 ) with respect to exchanging u and u is inherited from ( 2.22|) . 



B. Naive Perturbation Theory for the velocity field 

It is natural to seek solutions of Eq. ([2.17|) in the form 



u(r,t) = £u n , u n ocr n . (2.23) 



n=0 



Substituting this in the Navier-Stokes equation in the form (|2.17 ), and collecting terms with 



the same power in T we have a set of recurrence relations starting with 

m = * r-° > ( 2 - 24 ) 
u 2 = \g° * r-j + 1g° * rz = g° * r% . (2.25) 

The last equality follows from the symmetry of T which was discussed in the previous 
section. In fact, we are making full use here of the fact that we study a classical system. In 
the corresponding quantum mechanical theory the two terms in Eq. fl2.25| ) are distinct due 
to the non-commutativity of ui and u which become operators. 
The next orders have the form 

u 3 = G *r:- + lG *r^, (2.26) 
u 4 = g° * r;^ + g° * r »i , (2.27) 



9 



etc. In general, 



n/2-X 

u n = ]T G *^-^ 1 (neven), (2.28) 

m=0 
(n-3)/2 

u n = J2 G°*r;z- m+1 

1 

— I 

2 



+-G°*r:^:^ (nodd). (2.29) 



At this point one may substitute lower order u m in ( p. 28] ) and ( |2.29| ) to achieve eventually 
a solution in terms of uo and G°. For example 112 can be written as 

u 2 = G° * rii 7 ^^ (2.30) 
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FIGURES 



U ~ AAA i U o ~ AAA , 

* r: ~ -£ 

F ~ AAAAAA, F Q ~ AAAAAA, 

G ~ An n G Q ~ AATl . 



FIG. 1. Graphical notation for the perturabtion expansion. The symbol used are the following: 
short wavy lines stand for the fluid velocity (think about waves). A thin wavy line stands for uo, 
whereas a bold wavy line represents the full solution u(r,i). A straight line stands for the field 
p(r, t) that is only introduced in Chapter III. The Green's function, which is the response in the 
velocity to some force is made of a short wavy line and a short straight line representing the force. 
Again this symbol appears in thin and bold variants. The former stands for the bare Green's 
function (2.11), and the bold for the dressed Green's function (3.5). The vertex (2.18) is a fat dot 
with three tails. One straight tail belongs to the Green's function, and two wavy tails stand for 
velocities. A long wavy line will represent correlation functions of velocities. Again, the thin and 
bold variants are bare and dressed correlators respectively. 

For higher values of n this substitution becomes increasingly more cumbersome in ana- 
lytic form, and it is very helpful to represent it graphically. The notation used is shown in 
Fig. |l|. In Fig. |2|a we display the diagrammatic representation of the equation of motion, and 
in Fig. |^b and Fig. ^|c we show the diagrammatic form of Eqs. ( p. 24 ) and Q2.3C ). Examining 



the other graphs in Fig. |2] we see the great advantage of the diagrammatic representation. 
It is very easy to write the n'th order diagrams without going through the cumbersome 
analytic substitutions. The rules for drawing the diagrams are very simple: for the n'th 
order term we draw all the topologically distinct binary trees with n vertices, such that all 
the trunks are made of Green's functions and all the end branches are made of Uo's. Every 
portion of the tree that continues in a symmetric fashion gets a factor of 1/2. That is all. 
The existence of these simple rules for finding the n'th order term, which are independent of 
n, reflects the internal sturcture of this perturabtion theory, and it will allow us to perform 
resummations in a straightforward way. 

It should be understood that the perturbation theory described in this section is not 
expected to converge. Consider for example the ratio of U! over u . "Dividing" Eq. ( [2.24| ) 
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by mq leaves us with the order of magnitude estimate 



ui/uo ~ G° * r; uo ~ Vu /zA7 2 ~ u L/v = Re . 



(2.31) 



We see that we are effectively expanding in Reynolds number, and for high Re this is a 
terribly divergent series. The main problem in this expansion is the appearance of the 
molecular viscosity v in the denominator in (|2.31| ). It is expected that in a turbulent fluid 
the viscosity is strongly renormalized because of the interaction, yielding the so called "eddy- 
viscosity". We need to reformulate the theory such that the eddy- viscosity appears in the 
effective expansion parameter to make it of order 1. This is achieved in Sect.|2]. 

(a) 



u = 



+ 



(b) 


(c) 

2 3 


Uj = 

\ 1/2 

I 'a 


i V 

u 2 = \J " 2 




(d) 

3 4 


12 3 4 




i/r \/ 7/2 + 




U 3 = 1/2 J 

I 3 A 


I 3B 
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FIG. 2. Graphic representation of the equation of motion [Panel (a)] and the naive perturba- 
tion theory. The diagrams in Panels (b)-(d) are all the contributions appearing in ui — 114. Note 
that all the topologically distinct trees appear. These trees have a trunk that begins with the wavy 
line of a bare Green's function, and continues all the way via bare Green's functions. The branches 
are made of wavy lines of uo. The bifurcations carry a vertex T. The factors of 1/2 appear in all 
the bifurcations that have symmetry above them. The total factor in front of a tree is the product 
of these numbers. Every tree has been denoted with a letter and a superscript. The latter stands 
for the order in T n . For future purposes we have numbered the branches of the trees 2 A, 3 A, 3 B 
and 4 A 



C. Naive perturbation theory for statistical quantities 



1. Statistics 



As explained in the Introduction, our aim is to formulate a theory for ergodic averages. 



In the naive perturbation theory we can compute such averages by using the expansion (|2.23|) 
for u(r, t), and averaging products of uo(r,t) at different points and times with respect to 
different realizations of the random forcing f . In doing so we shall make full use of the fact 
that f is assumed to be Gaussian. On the one hand, this leads to a great simplification of the 
naive perturbation theory. On the other hand it seems theoretically dangerous. The physics 
that we are trying to describe is surely non-Gaussian. Are we not throwing away important 
aspects of the physics by using from the start a Gaussian forcing? The answer will be NO. 
We shall see later that the structure of the theory in its final resummed formualtion does 
not depend on this way of constructing the naive expansion. The physics is dominated by 
strong interactions via the nonlinearity, and a Gaussian forcing is sufficient to excite all the 
possible responses which the dynamics amplifies and redistributes to allow us to observe the 
full statistical structure. In fact, at the end of the calculation we are going to explore the 
limit in which the random forcing becomes vanishingly small. We shall show that this limit 
exists, and that in that limit the statistics of velocity fluctuations becomes independent of 
the statistics of the random force. Notwithstanding, it should be said here that it is not 
at all obvious that the procedure followed here is the most efficient one. It is possible that 
one can get a better converging theory by forcing the system with a force that is closer in 
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character to the statistics of the dressed system. We do not know how to do that, and this 
is one of the open problems that needs to be studied in the future. We shall return to this 
issue in Sec. 3 after the formal apparatus needed for its full appreciation is already at hand. 

We will denote averages with respect to the realization of the random force with single 
pointed brackets, in contrast with the double pointed brackets of Eq. (|1.1|). The Gaussianity 
of the random force means that 



(!(*)> = o, 

({(x 1 )~f(x 2 )) = D 12 , 
(f(xi)f (x 2 )f (x 3 )> = 

(f (Xx)f (X 2 )f (X 3 )f (X 4 )> = ^12^34 + ^13^24 + ^13^23 , 



^2.32) 



etc. All higher order correlations with odd number f vanish, and the correlations involving 
In factors of f have (2n — 1)!! contirbutions corresponding to all the possible pairings of f. 



(a) 



(h) 
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FIG. 3. The diagrams representing the mean velocity. Panel (a) represents the first order 
contribution in T, and is obtained from the simple gluing of the two branches of the tree 1 A in 
Fig. 2b. It carries the factor of 1/2 which originates from the symmetry of 1 A. Clearly, averaging 
either 2 A of Fig. 2c or any tree in Fig 2e results in a zero contribution due to the odd number of 
branches which carry a random force: (112) = (114) = 0. Panel (b) represents the result of averaging 
of the trees 3 A and 3 B in Fig. 2d. The diagran 3 Ai originates from gluing of the branches 1 with 2 
and 3 with 4 in 3 A. Diagram 3 yl2 comes from gluing either 2 with 3 and 1 with 4 or 1 with 3 and 
2 with 4 in 3 A. Correspondngly we gain a factor of 2 with 3 'A2 leading to a coefficient 1/4 instead 
of 1/8 in 3 A. The tree 3 Bi originates from 5 B after gluing 1 with 3 and 2 with 4 or 1 with 4 and 
2 with 3. A factor of 2 is gained. Lastly, diagram 3 I?2 results from gluing 1 with 2 and 3 with 
4 in 3 B. There is only one way of doing it, and the factor 1/2 remains. The general rule for the 
overall factor in front of a diagram is obtained as follows: count the number of vertices such that 
that exchanging the two branches emanating from them leaves the diagram invariant. Denote the 
numebr of distinct pairs of such branches by N. The overall factor in front of the diagram is 1/2 N . 



2. The mean velocity 

Of the sought statistical quantities, the easiest to obtain is the mean velocity, averaged 
over all the possible realizations of the random force. Since this random force is Gaussian, 
we have well denned statistics for the averaging process. We can apply the rules (|2.32|) to 
the average of the diagrammatic representation of u(r, t). We pair the Uo branches in all the 
possible ways, and glue the ends together, see Fig. || Every diagram with an odd number 
of u branches gives no contribution. Every diagram with 2n u branches gives (2n- 1)!! 
contributions which are obtained from all the possible binary pairings of random forces. The 
process is shown in Fig. [| The diagrams n A m and n B m for (u n ) in Fig. ^| result from the 
diagrams n A and n B for u n in Fig. |^. Note that in systems which are homogeneous and 
isotropic the mean velocity vanishes. Consequently, the sum of all the diagrams obtained 
in this fashion has to vanish. This will be used in our later developments. In a turbulent 
system with a space dependent mean velocity profile this set of diagrams will not vanish, and 
it will contribute also in other statitsitcal averages that we consider below. These diagrams 
will describe the interaction of the mean profile with the velocity fluctuations. 
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3. The Green's function 



Next we discuss the Green's function, which is the response of the velocity field to an 
external perturbation. The Green's function is defined as 

G a p(x,x') = i^5u a (x)/5ff 3 (x')' K j (2.33) 

where the notation $(■)/$(■) stands for the functional derivative. The meaning of this func- 
tional derivative is the following: solve the Navier-Stokes equations once with a forcing f 
and once with a forcing f + e5(x — x'). Then take the ratio [u(x, e) — u(x, 0)]/e in the limit 
e — > 0, and average over the ralizations of the random force. The principle of causality 
means that G(x,x') is zero for t' < t. This property will be used a lot in the sequel. 

The calculation of the functional derivative using the diagrams in Fig. ^ is straightfor- 
ward. Every diagram having n branches of uo contains a product of n f 's. The calculation 
of the derivative with respect to 8f(x') means via the chain rule that we get n contributions 
to 5u(x)/5f (x f ) . Every such contribution is obtained by dropping one of the branches of u , 
and replacing it by branch of G°. An example of how this procedure is done for diagrams 
2 A and 4 A in Fig. |2] is shown in Fig. f|a and Fig. respectively. 

The second step is obtained by averaging the diagrams for 5u(x)/5f(x') over realizations 
of f(x'). We can apply the rules Q2.32 ) to the averages of the diagram of 5u(x)/5f(x') in a 



graphical sense. This procedure is done in Fig. |5|. Fig. |5]a contains the two contributions to 
G 2 , and Figs. f|b and ||c show all the 19 diagrams for G 4 . The diagram denoted n A mp are 
obtained from the diagrams n A m in Fig. |4|. The diagrams denoted n B mp and n C mp originate 
from diagrams n B and n C in Fig. |^. Notice that every diagram has an "entry" which is the 
root of the tree in Fig. and an "exit" which is the branch G° which replaced a branch of 
u . There is a unique path between entry and exit which is composed of a chain of G° 's. 
We shall call this path the "principal path" of the diagram. Notice that the entry always 
begins with a wavy line, whereas the exit ends with a straight line. In fact, all the diagrams 
appearing here can be drawn without reference to the explicit derivation described here 
using simple topological rules. However, since the diagrams described here are not in their 
final form we defer the discussion of the appropriate rules for a later moment. 
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FIG. 4. Typical contribution to the diagrammatic representation of 5u(x)/5i which originates 
from the trees 2 A and 4 A in Fig. 2c and 2e. In panel (a) we show the diagrams originating from 2 A 
and in panel (b) those originating from 4 A The diagram 2 A2 is obtained by differentiating with 
respect to f(x') at position 2 or 3 in 2 A. A factor of 2 is gained. The diagram 2 A\ comes from 
differentiating with respct to f(x') at position 1. The diagram A± is obtained by differentiating 
with respect to f(x') at position 4 or 5 in 4 A. A factor of 2 is gained. The diagrams A A^ , A A2 and 
4 ^4i come from differentiating with respect to f(x') at positions 3, 2 and 1 respectively. Again the 
factor 1/2 remains at vertices that have symmetry above them. 
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(c) 




FIG. 5. Diagrams for the Green's function. Panel (a): The contributions to G2 originate 
from the trees shown in Fig. 4a by gluing two uo branches in the only possible way. Panel (b): 
The contributions to G4 originating from the trees A\ — A± shown in Fig. 4b. The diagrams 
4 ^4i — 4 A43 come from 4 ^4 by gluing 1-3 and 2 - 4; 1 - 4 and 2 - 3 ; 1 - 2 and 3-4 repsectively. 
The diagrams 4 Aji come from 4 A3 by gluing 1-4 and 2 - 5 or 1 - 5 and 2 - 4. A factor of 2 is 
gained. 4 A%2 is obtained by gluing 4 A3 at positions 1-2 and 4-5. The diagram 4 ^2i comes from 
4 ^2 by gluing 1-3 and 4-5. 4 ^22 is obtained by gluing 1-4 and 3-5 or 1-5 and 3-4, gaining a factor 
of 2. Finally 4 in comes from 4 A\ by gluing 2-4 and 3-5 or 2-5 and 3-4. A\i obtains from 4 A\ by 
gluing 2-3 and 4-5. (c) Diagrams for G4 originating from 4 B and 4 C in Fig. 2e. Note that again 
all the numerical factors follow the same rules as observed in Fig. 4. 
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4- The 2-point velocity correlation function 



The 2-point velocity correlation function F(x,x') is defined as 

F aP {x,x') = {u a {x)u p {x')) . (2.34) 

The calculation of this quantity is again based on the diagrammatic expansion shown in 
Fig. 0. To obtain F m+P of order Y m+P we need first to take a contribution u m and a 
contribution u p and mutiply them together. The n-th order F n is obtained as a sum of 
all F m+P such that m + p = n. In the second step we average over the realizations of f . 
As before all odd order contributions vanish because they contain an odd power of f . For 
example, 

F (x,x') = F 0+0 (x,x'), 

Fx(x, x') = F 1+0 (x, x') + F 0+1 (x, x') = 0, (2.35) 
F 2 (x, x) = F 2+0 (x, x') + F 1+1 (x, x') + F 2+0 (x, x') , 



where 



F m+P (x,x') = (u rn (x)u p (x')) . (2.36) 



To obtain the diagrammatic representation of this procedure we take one tree for u m (x) 
and one tree for u p (x') from Fig. 0, and average the product of these trees according to the 



rules ( 2.32 ). In other words we need to pair the branches of u in all the possible ways, and 
to glue them as discussed before in computing (u) and the Green's function. The procedure 
is shown in Fig. |6|. Note that every diagran has a uniquely defined "principal cross section" 
which arises from the gluing of the two trees u m (x) and u p (x). We denote it in the diagrams 
as a vertical broken line, and we draw the reader's attention to the fact that the principal 
cross section cuts through correlators only, and not through Green's functions. This fact is 
used later in the process of resummation. 

This is the end of the naive perturbation theory. As we saw, it is an expansion in powers 
of Re, and we are going to partly resum this expansion to develop a reformulation with a 
better expansion parameter. We reiterate here that the procedure at this point seems very 
dependent on the properties of the noise, since we used the rules (|2.32|) time and again. We 



defer further discussion of this issue until after the resummations of various sorts, when we 
can take the limit f — > with impunity. 



19 



3. RESUMMATIONS 



In this chapter we discuss the resummation of the naive perturbation theory that was 
developed in Chapter [2|. We begin with the mean velocity and its role in the resummation 
of the Green's function and the correlator. 
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FIG. 6. Diagrammatic representation of the 2-point velocity correlation function. Panel (a): 
Fi+i is obtained from the the trees Ui(x) and ui(x'). Diagram 1 is obtained by gluing 1-2 and 
3-4, and diagram 2 by gluing 1-3 and 2-4 or 1-4 and 2-3, gaining a factor of 2. Panel (b): ^0+2 is 
obtained from the trees for u(x) and 112(0;'). Diagram 1 is obtained by gluing 1-2 and 3-4, while 
the diagram 2 comes from gluing 1-3 and 2-4 or 1-4 and 2-3, gaining a factor of 2. Panel (c): Some 
typical diagrams contributing to F4. Note that all the diagrams have a "principal cross section" 
that separates trunks belonging to the left and to the right trees. This cross section runs through 
2-point correlators, never through Green's functions. In contrast with the diagrams for the Green's 
function that have wavy entries and straight exits, here we have wavy entries on both sides. The 
rules for the numerical factors of the diagrams remain the same. 

A. The resummation of the mean velocity 

In Section |2|B we discussed the diagrammatic series for (u) (see Fig. |3|) , and commented 
that it has to sum up to zero when (u) = 0. In discussing the diagrams for the Green's 
fucntion Fig. [| and for the correlator Fig. |B| we encounter again the same type of diagrams 
that appear in (u). For example in Fig. [5] the diagrams 2 A U , 4 t4 32 , a A 2 i, 4 v4 n , 4 B 2 i, 4 B 31 , 
and 4 Cn all have a fragment which is the diagram for (ui) in Fig. ^a. In addition, the 
diagram 4 Au has a fragment which is identical to 3 Bi in Fig. |3|. The diagram 4 Au has a 
fragment like 3 B 2 . The diagram 4 C 2 i has a fragment like 3 A 2 , and lastly 4 C*22 exhibits a 
fragment like 3 A\. 

All these diagrams that contain fragments belonging to (u) have a common feature. To 
see this denote the part of any diagram that contains the principal path between entry and 
exit as the "body" of the diagram. Any fragment that can be disconnected from the body 
by cutting one Green's function is called a "weakly linked" fragment. All the diagrams that 
we discussed in the previous paragraph have weakly linked fragments. Consider now all the 
diagrams in Fig. |5| that have two Green's functions in their principal path. These are the 
diagrams 2 An, A An, A Ai 2 , A C 2 \ and 4 C*22. The sum of the weakly linked fragments of these 
diagrams is exactly (ui) + (113) as can be seen in Fig. || If we consider all the higher order 
diagrams for G which have two Green's functions in the principal path, we find that their 
weakly connected fragments furnish all the remaining diagrams in the series for (u). The 
coefficients in front of all these fragments is the same as the coefficient in the series for (u) 
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since it is determined uniquely by the local topology of the fragment, independently of the 
position that the fragment occupies in the mother diagram. Accordingly all these diagrams 
with two Green's functions in the principal path sum up to zero. The same story repeats 
for all the diagrams that contain a weakly linked fragment. For example the diagram 4 ^2i 
in Fig |5|b is the first in the series that exhibits a weakly linked fragment that eventually 
will be resummed together with diagrams that have the same body, but higher order weakly 
linked contributions that sum up to zero. The general conclusion is that all the diagrams 
that have at least one weakly linked fragment sum up to zero and need not be considered 
further in the resummed theory. 

Next we discuss the appearance of (u) in the series for the 2-point correlation fucntion. 
In this series we find a new type of diagram, like diagram 1 in Fig. ||a. These are unlinked 
diagrams which are obtained from averaging the left and the right tree separately. Such 
diagrams contain no correlators that cross the principal cross section. Obviously such dia- 
grams will resum to (u) 2 which is zero. In addition we have diagrams with weakly connected 
fragments. In the context of the diagrams for the 2-point correlator we define the body of 
the diagram as the part that contains the two entries (the roots of the original trees), which 
are denoted by "x" and "x"' in all the diagrams in Fig. ^. A weakly linked fragment is a 
fragment that can be disconnected from the body of the diagram by cutting off one Green's 
function. An example of such a diagram is diagram 1 in Fig. |6|b. As in the case of the 
Green's function, the infinite sets of weakly linked fragments with the same body resums to 
(u) = 0. From this point on we therefore discard all the diagrams that have at least one 
weakly connected fragment. 



B. The Dyson resummation for the Green's function 

In this section we discuss the Dyson line resummation of the series for the Green's 
function. To this aim we classify the diagrams into three classes. The first class consists of 
diagrams with weakly linked fragments which are all discarded. The other two classes are 
designated as follows: 

I. Principal path reducible diagrams. These are diagrams that can be split into two 
disjoint pieces which contain more than one Green's function, by cutting one bare Green's 
function that belongs to the principal path. An example of such a diagram is 4 ^43 in Fig. ^b. 
This diagram fall into two parts by cutting Gq(x2,xs). 

II. Principal path irreducible diagrams. These are the diagrams that cannot be split as 
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described in I. All the other diagrams in Fig. [5] that do not have weakly linked fragments 
are principal path irreducible. 

All the principal path irreducible diagrams, except G (x,x') itself, share the property 
that they start with a bare Go(x,Xx), they end up with a bare Go(x 2 , x'), and in between 
they have a principal path irreducible structure, say S(xi,x 2 ). The sum of all these ir- 
reducible structures is defined as the XI operator, which will be shown to contain all the 
information about the turbulent eddy viscosity. Using the diagrams in Fig. [| we develop the 
diagrammatic expansion for S(x, X\) which is shown in Fig. [7[ With the help of H,(xi,x 2 ) 
we can say that the sum of all the principal path irreducible diagrams is 

Sum of all irreducibles (3.1) 
= G (x, x') + G (x, xi) * E(xi, x 2 ) ■ G (x 2 , x') 



where the star "*" and dot "•" products are as defined in the vertex Eq. (|2.17Q . The operator 



S(xi,X2) starts and ends with a vertex, and it therefore connects with "*" to the preceding 
G (x, Xi), and with a "•" to the following G (x, Xi). Further summation of this series will 
be discussed soon. 

The principal path reducible diagrams, of which we have only one representative in 
Fig. |5|, can be also resummed using the operator 22). Note that the reducible diagram 

4 A« has a structure such that that between x\ and x 2 there exists the first contribution to 
~S(xi,x 2 ), and after the point x 2 we see a fragment that is identical to 2 A 2 i, which is the 
first nonlinear contribution to the Green's function. The main observation, which is due to 
Dyson, is that higher order reducible diagrams will have between x\ and x 2 all the higher 
order terms in x 2 ), and than after x 2 we will have all the other nonlinear contributions 
to the Green's function G(x 2 ,x'). Therefore, we can write 

Sum of all reducibles (3.2) 
= G (x, x x ) * S(xi, x 2 ) ■ [G(x 2 , - G {x 2 , x')} 

Again we made use of the fact that all topologically possible diagrams appear in the series, 
and that the numerical weight of each fragment is only determined by its local symmetry, 
independent of its position in the mother diagram. 
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FIG. 7. Panel (a): Diagrammatic series for the £(xi,x 2 ) operator. All the diagrams to order 
r 4 are shown. Every diagram is denoted by the same notation of the diagram it derives from in 
Fig. 5. Panel (b): Diagrammatic series for the <$>(xi,X2) operator. All the diagrams to order T 4 
are shown. The numbering of the diagrams is referred to in the section on line resummation. 



Adding together ( |3.1| ) and ( |3.2| ) we get the Dyson equation for the Green's function 
G(x, x) = Gq(x, x 1 ) + Gq(x, x\) * S(xi, x 2 ) • G(x 2 , x) (3.3) 



The formal solution of this equation in operator form is 



G = [l-Go*S]- 1 G t 



(3.4) 



or, more explicitly 



G=[G -S]- 1 



iP 



d/dt - z/V 2 - PS 
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(3.5) 



where we insert the definition of Go from the Eq. [2.11 ) has been used. We see now that the 
S operator serves as a renormalization of the viscosity. The full appreciation of this fact 
will become clear in subsection D. 



C. The Wyld resummation for the 2-point correlation function 

In this section we discuss the Wyld line resummation of the series of the 2-point correla- 
tion function. We need again to classify the diagrams in order to see clearly how to proceed. 
Firstly, all the diagrams that are unlinked or have weakly linked fragments can be discarded. 
Next we consider all the diagrams that have only one wavy line that crosses the principal 
cross section. The diagram 2 in Fig. [5]b and the diagram 3 in Fig. |6]c are such diagrams. In 
Fig. [8] we show a typical 8 th order diagram of this type. 




(b) 



AAA/1 — 




FIG. 8. Panel (a). One of the diagrams that have one wavy line at the principal cross section 
which appears in the series for the 2-point correlation function. Panel (b). The representation of 
the lowest order contribution to the 2-point correlator 

Recall now that the wavy line, which is a correlation function of two u 's, {u (x)u (x')) , 
can be interpreted as Gq(x, x\) * (f(xi)f (xz)) * Go(x2,x') (see Fig. |8|b). Thus, on the two 
sides of the (f(xi)f(x2)) we have two typical contributions to the diagrammatic series of 
G that was discussed in the last section. On both sides of the principal cross section we 
get, excluding (f(ari)f (a^)), portions of diagrams that have an entry and an exit like all the 
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diagrams in Fig. |5|. Thus the diagrams 2 in Fig. |6|b and 3 in Fig. ^|c have fragments which 
are identical to 2 A 2 \ in Fig. ||a. In the example of Fig. [8] the two fragments are identical with 
diagrams 4 -Bi2 and 4 A 43 in Fig. |5|. All the other diagrams with only one wavy line across 
the principal cross sections will have all the other contributions to the series of G on both 
sides of the (f {x-\)f (x 2 )) correlation. Thus we conclude that 

Sum of all diagrams with one 
wavy line at the cross section 
= G(x, xi) * (f (xi)f (x 2 )) * G(x 2 , x') . (3.6) 

Next we consider the diagrams with two or more wavy lines crossing the principal cross 
section. Examples are diagram 2 in Fig. |6|a, and diagrams 1,2,4, 5 and 6 in Fig. ^|c. Of these 
diagrams consider first all those that cannot be split into two parts (containing more than 
a single Green's function) by cutting off one Green's function. All these diagrams have two 
entries via straight lines on the right and the left, and a structure, say 0(xi, x 2 ) in between. 
The sum of all these 4>(xi,x 2 ) contributions is defined as the Q(xi,x 2 ) operator. In Fig. [7|b 
all the contributions of $(xi, x 2 ) up to 4 th order are exhibited. 

All the rest of the diagrams that have two or more wavy lines crossing the principal 
cross section are analyzed as follows: start from the principal cross section, and move left 
until you meet for the first time a Green's function that can be cut such that the diagram 
splits into two disjoint parts. Denote the coordinate of the vertex at the exit of this Green's 
function as x\. Start again from the the principal cross section and move the right until 
you meet another Green's function that can be cut such that the diagram splits into two 
parts. Denote the coordinate at the vertex at the exit by x 2 . This notation is shown in all 
the relevant diagrams of Fig. |6|. The fragment to the left of X\ and to the right of x 2 is a 
typical fragment that appears in the series for G. For example in diagram 2 Fig. |G|a we have 
a Go(a;,:Ei) on the left of X\ and Go(x2,a;') on the right of x 2 . In diagram 1 in Fig. |]c we 
have a fragment 2 A 2 i (Fig. |J) to the left of X\ and again Go to right of x 2 . In the diagram 
2 in Fig. ||c we have the inverted situation. Obviously, in the higher order diagrams we will 
find all the remaining contributions to G. We thus conclude that all these diagrams which 
have two or more wavy lines crossing the principal crosss section may be resummed to 

Sum of all diagrams with 2 and more 

wavy lines crossing the cross section (3-7) 
= G(x, xi) * x 2 ) * G(x 2 , x) . 
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Summing up Eqs. (|3.6| ) and ( |3.7| ) we recover the Wyld equation 



F(x,x') (3.8) 
= G(x, xi) * [(f (xi)f (x 2 )) + x 2 )] * G(a? 2 , a/) . 

This equation, after the resummation of the series for <&, will allow us to discuss the limit 
f->0. 



D. Line resummation 

Up to now our perturbation series is in orders of the Re number, as explained in Sect|| 
The Green's function and the 2-point correlators which appear in the diagrams for $ and 
S are the bare ones. In this section we perform the so called "line resummation" which 
results in diagrams in which all the propagators are dressed ones, and at the same time the 
perturbation parameter becomes of the order of 1. 

To accomplish this resummation we introduce a few notions. First is the notion of the 
"end points" of a diagram. The end points are defined as the entry and exit in the case of the 
dagrams for X, and the two exits in the case of diagrams for <fr. Next we define "one eddy 
reducible" diagrams and "one-eddy irreducible" diagrams. These types of diagrams play 
a topologically similar role to the "one-particle reducible and irreducible" diagrams in the 
standard Feynman diagrams. The definition is as follows: Name as a "one-eddy reducible" 
fragment any fragment that can be completely disjoined from the part of a diagram that 
contains the end points by cutting any two lines, be whether Green's fucntions or correlators. 
The part that contains the two end points will be called the bulk of the diagram, which may 
or may not be the same as the body of the diagram as defined in SectiorJ3~M The convention 



will be that the fragment contains the two lines that have been cut. In other words, we cut 
lines near the vertex that belongs to the rest of the diagram. It follows from this convention 
that one-eddy reducible diagrams are at least of order T 2 . Define now 

I. One-eddy reducible diagrams: diagrams that contain at least one one-eddy reducible 
fragment. 

II. One-eddy irreducible diagrams: diagrams that contain no one-eddy reducible frag- 
ment. 

Examples for one-eddy reducible diagrams are diag rams ^^42, ^^22? ^-^22 and ^B-$2 in 
Fig. [7|a for X, and diagrams 5, 6 and 7 in Fig. |7|b for The rest of the diagrams in Fig. [7] 
are one-eddy irreducible. 
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We can classify one-eddy reducible fragments into two classes, called as class F and class 
G. Class F has two wavy ends (like the fragment of A A 22 ) and class G has one wavy and 
one straight end (like the fragment of 4 A 42 ). An observation of principal significance is 
that the sum of all the possible fragments of class F is precisely the nonlinear correction 
to the 2-point correlator F, whereas the sum of all the possible fragments of class G is the 
nonlinear correction to the Green's function G. As an example the fragment of 2 A^ 2 is the 
first nonlinear correction to G which is 2 A 2 \ in Fig. |5|a. The fragments of 4 B 22 in the series 
for X, Fig. ^|a and the fragment of diagram 5 in the series for Fig. [?|b are the same, and 
are the nonlinear correction 2 for F in Fig. ^|a. The fragment of the diagram 4 A 22 is the 
diagram 2 in Fig. ^p. 

Consider now the irreducible diagram 2 A 2 \ in Fig. [7]a. By summing it together with the 
one-eddy reducible diagrams whose bulk part contains either one bare correlator or one bare 
Green's function we obtain the diagram 1 in Fig. |S]a. Similarly, every irreducible diagram 
in the series for S and $ acts as the beginning of an infinite series of one-eddy reducible 
diagrams that renormalize all bare propagators to dressed ones. In Fig. ^] we show all the 
contributions up to 4'th order of the resulting series. This (very important) step in the 
theory is available due to two facts: the first is that all the topologically allowed diagrams 
are present, and the second is that the weight of each fragment is determined locally by its 
symmetry properties. 

We note that at this point one can state the following rules of how to draw any diagram 
in the renormalized series for X and 

I. Only diagrams with an even number of vertices appear. 

II. The series for S and <fr contains all the topologically allowed diagrams made of 2n 
vertices with one straight and two wavy legs. The difference is that the diagrams for S must 
have an entry with wavy line and an exit with a straight line, whereas diagrams for <fr have 
two exits with straight lines. 

III. One needs to exclude all the one-eddy reducible diagrams, and also all diagrams 
that include closed loops of Green's functions. The latter are zero by causality. They never 
appear in the method of derivation discussed above, but they do appear in the path integral 
formulation that is discussed in Section [|. They vanish also there, of course. From these 
properties it follows that the diagrams for I] have a unique principal path, and that the 
diagrams for $ have a principal cross section. 

IV. The coefficient in front of a diagram is determined by the symmetry properties. 
Every flip of legs that leaves the diagram invariant contributes a factor of 1/2. However, the 
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one-eddy irreducible diagrams for I] cannot have any element of symmetry, and therefore 
all their weights are unity. The irreducible diagrams for F can have at most one element of 
symmetry, and they will have a weight of 1/2. Examples are diagrams 1,5 and 6 in Fig. [|b. 

We can show now that the effective coupling constant is of the order of unity. First we 
discuss how many correlators, propagators and vertices appear in a 2n'th order diagram for 
X and In a symbolic fashion we can write that 

S 2 „ ~ H 2 K n ~ l , S 2 ~ (T 2 GF) ■ (3.9) 

$ 2 „~$ 2 A n ~\ $ 2 ~r 2 F 2 (3.10) 

where the new expansion parameter A is of the order of 

A ~ T 2 G 2 F . (3.11) 

To estimate the order of magnitude of A we notice that G, according to the Dyson equation 
fl3.4|) , is of the order of S -1 . Using this in ( |3.11| ) we get 

a.E!^^. ,, 12) 

where Eq. ( |3.9|) has been used. We see that the line-renormalization has shifted our expan- 
sion parameter from the order of Re to the order of 1. 

At this point one usually considers the possibility of vertex renormalization. In our 
context this can be achieved by looking at the 2-eddy reducible fragments that can be cut 
from the bulk by severing three legs. It will turn out that our perturbation series has the 
deeply non-trivial property that that this step is not needed in the final presentation of the 



theory. This important issue is discussed further in [21 



We thus end with the Wyld equation which can be rewritten as 

F(x, x') = G(x, xi) * [(f(xi)f(x 2 )) (3.13) 
+ {cf){x 1 )cj)(x 2 )) + §(x 1 ,x 2 )] * G(x 2 ,x') 

where we have made use of the fact that the thermal noise f is uncorrelated with the stirring 
force 4>. We are going to seek a solution of ( |3.13 ) for which the energy of the turbulent 
motion per mode is much larger than the order of k B T. For such solutions the effect of 
the thermal noise on the observed statisitcs of the velocity fluctuations will be negligible, 
and we can omit (f (xi)f (x 2 )) from (|3.13|) with impunity. On the other hand, the role of the 
stirring force (4>(xi)<fi(x 2 )) can be very important. In general the observed statistics depends 
on how the fluid is stirred. The choice of stirring that models the boundary condition of 
experimental turbuelence is discussed in [|16|] after performing some of the analysis on the 
operators S and <3>. 
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E. Intuitive meaning of the Dyson and Wyld equations 



The physical significance of the Dyson equation (|3.3|) is that the dressed response, which 
is a function of r, r' and t — t', is determined by hydrodynamic interactions involving inter- 
mediate points. For example, the response to forcing at r' has one direct contribution at 
r, which is Gq(x, x'). However, for any finite time the response to forcing at r' is mediated 
by interactions at points ri which via S appear at point r. In the one loop approximation, 
which is diagram 1 in Fig. [| I] itself has a Green's function that mediates directly between 
ri and r. In higher order contributions to S there are sequential contributions due to forc- 
ing at r' that are mediated by responses at r 1; r 2 ... until r is reached. X represents the 
dressed response which is the sum of all these sequential responses at multiple intermediate 
sets of points. The Green's functions that mediate intermediate points are weighted by the 
correlators of velocity differences between these points; if these correlators are small, the 
contribution of the Green's fucntion to the total response is also small. 



(a) 




FIG. 9. Diagramatic representation of of the operators X and <3? in the line-renormalized theory. 
The symbols for the Green's function and the 2-point correlators are bold; they are renormalized 
to all orders. 
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The intuitive understanding of the Wyld equation is also straightforward. From the 
equation of motion ( 2.10|) , written schematically as 



u =wra¥- |u ' v|u+f} ' (314) 

we see that the nonlinear term [u ■ V] u can be understood as an "additional noise" in the 
equation of motion. It is natural to expect that the double correlator of the forcing will have 
a nonlinear contribution of the type 

(fWOf^))^ (3.15) 
~ ([u(ri,*i) • Vi]u(ri,*i)[u(r 2 ,* 2 ) • V 2 u](r 2 ,t 2 )) . 

Indeed, the mass operator $ Q/ a(r, r', t — t') can be written exactly as 

$ a/3 (r,r',t-t') (3.16) 
= q t g r , F ajl 36(r, r, r', r', t, t, t', f) , 

where F aif3 s is the fourth order correlation of u. Thus Eq. (|3~8|) can be understood as the 
result of squaring Eq. (|3.14j) and averaging, up to the dressing of the Green's functions. The 



content of the Wyld diagrammatics is the dressing of the bare Green's function that appears 
in Eq . (|3.140 ) and the representation of the fourth order F ai/3 s in terms of second order F a p. 
The dressing of the Green's function appears from the cross terms between [u • V]u and f 
in ( |3.14j ). Note that the analytic form of the Q a/ 3 in the 1-loop order, diagram 1 in Fig. ^b, 
follows directly from the Gaussian deomposition of F a ^g. 
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<v, v 2 v,> = i 



+ ( /VWW( B ) 





<3= < 





FIG. 10. Diagrammatic representation of the 3-point velocity correlator. Panel (a): There are 
3 types of dressed vertices. Panel (b): The diagrammatic expansion for the three vertices appearing 
in (a). The expansion is in terms of dressed propagators and bare vertices, and we show all the 
diagrams up to fifth order in the bare vertex. 

F. 3-point and higer order velocity correlation functions 

To calculate the n-point velocity correlation function we need to take n trees from Fig. |2| 
and average their product. Diagrammatically it means that we pair branches in all the 
possible ways, and glue. Next we need to resum all the unlinked diagrams (that sum up to 
zero), and all the weakly linked diagrams (with the body defined as the structure having 
n entries) which also give zero contribution. Next we perform a line-renormalization. The 
procedure is identical to the one described above, and we leave the details to the interested 



reader. The final result for the 3-point correlation function is shown in Fig. |10[ There 
appear three types of vertices that we denote as A,B and C. Their diagrammatic series is 
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shown in Fig. |T0|b. The series in 10b is in terms of bare vertices, since, as we said before 
there is no need to renormalize the vertices in this theory (see ||21||). 



e=0 

<t>=o 

FIG. 11. The condensed representation of the operators E and <3? in terms of the dressed 
vertices defined in Fig. 10b. 

It should be noted that the objects A, B and C can be used to condense the series for 
the operators S and <fr. This is shown in Fig. [TTJ. The formal elegance of this condensed 
presentation is of no known use in the analysis. 



4. FUNCTIONAL INTEGRAL FORMULATION 



A. Introduction 



In this chapter we discuss the path integral formulation of the statistical theory of Navier- 
Stokes turbulence. The path integral formulation does not lead to different diagrammatics 
than the one described and discussed in the last chapter. Its main advantage is two-fold. 
On one hand it allows a very compact representation, and on the other hand it offers the 
ability to effect an infinite partial resummation by changing variables in the path integral. 
This latter property is the one that will turn out to be very useful for us, allowing us to 
explore new renormalizations of the diagrammatic series which are different from the Wyld 



and Dyson resummations, see [21 



B. Setting up the generating functional 

The starting point for the statistical description of turbulence in the functional approach 
is the same equation of motion used in the Wyld direct diagrammatic expansion: 
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du/dt + (u • V)u - uV 2 u - Vp = f . 



(4.1) 



with the same properties of the random force f . As before, we need to find the statistical 
quantities by averaging over all the realizations of the random force f . 
Introduce the Navier Stokes functional of u, J\f{u(x)} as 



Af{u(x)} = du/dt + (u • V)u - zA7 2 u - Vp . 



(4.2) 



Given any functional of u, say A4{u(x)}, we can calculate its value for a given realization 
of u which in turn is determined by some realization of the random force f . This is done 
using the formal solution of Eqs. (|4.1|) - (|4.2|) 



according to 



u(x) = AT 1 ^)} 



M{u(x)} = M{Af^{f(x)}} . 



The RHS of this equation can be rewritten in the equivalent form 

M{N'~ 1 {f(x)}} 

M{u{x)}5\u{x) -Af-^fix^Duix) . 



(4.3) 



(4.4) 



(4.5) 



Here the functional integration is interpreted as the continuum limit 



M 

Du(x) = lim TT / d 3 u(xi)5\V ■ u(xj)l 

m-^oo- LJ - J 

i=l 



(4.6) 



where {x{\ = {r 4 , t { } and i runs over a discrete space-time grid and M is the total number 
of points in the 4-dimensional grid. The delta function <5[V • u(xj)] limits the integration to 
divergenceless contributions. The delta functional in ( |4.5| ) is defined as the continuum limit 

(4.7) 

M 



u(x)-Af- 1 {f(x)} 

M 

iim n / slu^-Ar- 1 ^)} 



M^oo J , 
1=1 



Note that Eq. ( [4.5[ ) is but a cumbersome generalization to the continuum of the obvious 
equation 



f( x ) = J f(y)$( x - y) d y 



(4.8) 



Remembering that 
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8[g 1 (v) ~x] = (dg / dx)8[g(x) - y) 
one may generalize to the continuum limit according to 
8[u(x) -Af-^iix)}] = J{u}5[Af{u(x)} - f(x)], 
where the Jacobian J{u} is given by 

SM a {u(x)} 



J{u} = det 



8u p (y) 

On the grid this Jacobian is the determinant of the 3M x 3M matrix 

5Af Q {u(xi)} 



5up(xj) 



a/3 



| + <»-V)-„V> 



+ dpu a j8(xi - Xj) . 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



for Ai at a 



In the Appendix we demonstrate that this determinant is unity. Thus Eq. 
given realization of f can be rewritten as 

M{u(x)\f(x)} = J M{u'{x)}5[N{u'{x)} - f(x))Du'(x) . (4.13) 

Next we need to consider the average of M. with respect to the realizations of u, (which in 
turn are determined by the realizations of f). As in the Wyld approach, it will be assumed 
here that f has Gaussian statistics. The average of any functional ^{f} can be written as 
the Gaussian functional integral 

(4.14) 



(7>{f}> = - / Df(x)V{f} 



x exp 



1 



- / f a (x)Dj(x - y)f p (y)dxdy 



The covariance matrix D a p is determined by the condition that for the case V{f} 
fa(x)f a (y) Eq. (|4.14j) would lead to the equation 



(fa(x)fa(y)) = D aP (x - y) 



The partition sum Z\ is 



Zi= Df(x)exp \- - / f a (x)D a l(x - y)f p (y)dxdy 



(4.15) 



(4.16) 



In perfoming this average for the functional A / i{u(x)|f (x)} of ( |1.13| ) it is convenient to 
rewrite the delta functional first in an exponential form such that all the / factors appear 
in the exponential. We use the representation of the 1-dimensional delta function 
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5 ( z ) = I 7T exp(-ipz) 

J -oo AH 

at every point of the grid. Thus 

1 



(4.17) 



5 



Af{u{x)} - f{x 



lim 



exp 



M 



Pa \%i 



i=l 



M^oo (27r) M 
N a {u(Xi)} - f a (Xi) } 



M 

[Y[dPc 

i=l 



X; 



(4.18) 



We average now Eq.( |4.13| ) using the general Gaussian recipe ( |4.14j ), representing the 
delta function as the continuum limit of ( 4.171) . The result reads 



(M{u{x)}) = — / Df(x)Du'(x)DpM{u'(x)} 
x exp U J dxp(x) ■ [N{W(x)} - f(x)] J (4.19) 
x exp |-- J f a (x)D^(x - y)f p {y)dxdy)\ . 

We remind the reader that our path integrals are limited to divergenceless contributions, cf. 
( pj~6l) . In Eq. ( p^|) Z 2 is the partition sum which is the RHS of flQDp without M{u(x)}. 
Note that it serves to cancel the formally divergent factor {2n) M . Also note that p appears 
only in scalar products with divergenceless fields, and any projection onto vector fields that 
are not divergenceless cancels from Dp with the same Dp that appears in Z 2 . It will be 
convenient to take p to be divergenceless from the start, by defining Dp similarly to Du in 
Eq.flOD. 

We can perform now the Gaussian integration over Df . One way of doing it is to change 
the f variables using a similarity transformation that diagonalizes the matrix D, then to 
perform the Gaussian integral at each space-time point independently, and lastly to rotate 
back with the similarity transformation. The final result is 



(M{u(x)}) = i J Du'(x)DpM{u'(x)} 



(4.20) 



x exp 



ji / p a (x)Af a {u' (x)}dx 



~ J p a {x)D c 



It is quite evident that (|4.20|) furnishes a starting point for the calculation of any correlation 
function (u(a;i)u(x2)---u(x n )). We recall however that the theory calls also for the calculation 
of the response or the Green's function. We show now that this calculation is also available 
from Eq.( [4.20| ). To see this imagine that we add to the force f some external deterministic 
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component h(x) that makes (u(xi)) non-zero. The addition of this component changes 
N a {u(x)} in the exponent in (|4.20|) to N a {u(x)} — h(x). Consider now the response (|2.33 ) 
which in this case can be represented as 

.5{u a {x{)) 



G a /3{Xi,X2j 



5h /3 {x 2 ) 
For a finite h 



(4.21) 



2 

Using this we can compute (|4.21|): 



(«a(^i)) = \\ Du(x)Dpu a (x 1 ) exp {i J pp[Afp{u(x)} 
1 ( 

h p (x)]dx - - p a (x)D a/3 (x - y)p p (y)dxdy} . (4.22) 



G a p(xi,x 2 ) = — J Du(x)Dpu a {x 1 )p l3 {x 2 ) (4.23) 



x exp{i / pj3(x)[Aff3{u(x)}]dx 



- ^ J Pa{x)D aP (x - y)pp(y)dxdy} = (u a (xi)pp(x 2 )) 
In the same way one can compute non-linear response functions like 



h _^ = (uaix^ppix^p^xs)) 



G Nh (x u x 2 ,x 3 ) (4.24) 
= 5 2 (?/ a (xi)) 
8hp{x 2 )8h 1 {x 2 ) 

etc. It is obvious that every additional functional derivative will appear as another factor 
of p in the correlator. In general it is useful to introduce the generating functional 

Z(l, m) = ^exp J dx[u(x) ■ l(x) + p(x) ■ m(x)]^ . (4.25) 

All the needed statistical averages can be obtained as a functional derivative of Z(l, m), 
taken at 1 = m = 0, for example 

dl a {x 2 )dlp{x 2 ) 
5 3 Z(l m) 

(uaixjupix^u^xs)) = ' , (4.27) 

dl a {x 2 )dlp{x 2 )bl^{x 3 ) 

etc. Similarly one can compute any kind of Green's function. For example, 

G a ,(x u x 2 )= f Z(I ' m) (4.28) 
dl a (x 2 )dmf3{x 2 ) 

G nh {x u x 2 ,x 3 ) = 5la ( X2)5lp fa )5rrhM > ( 4 - 29 ) 



.c^M^OM 2 ^)) _ <5 3 Z(1, 



111 ) 



8hj(x s ) 5l a (x 2 )5lj3(x 2 )Sm 1 (x 3 ) 
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(4.30) 



etc. Finally we express the generating functional Z(l, m), with the help of Eq. (|4.20|) , in the 
form of a functional integral: 

Z(l,m) = — [ Du(i;)Z)p(x) (4.31) 

Zj J 



x exp jil + J dx[u(x) ■ l(x) + p(x) • m(x)]| . 



The quantity / is referred to as the effective action, and for future work it is useful to divide 
into two parts, the one quadratic and the other triadic in the field p and u: 

I = h + 4* , (4.32) 
dx\p a -^ - vp a \7 2 u a 

+ I p a (x)D afS (x - y)pp(y)dxdy, (4.33) 

r^-lf^M^ (4,4, 

1 f dq 1 dq 2 dq 3 



2 J 2tt 



,12 



-PaTap-yiq, qi, q2)up(q 2 )u J (q 3 ) 



The last line follows from the definition of the vertex in the q representation, see Eq. ( |2.22 ). 



Note that we did not display a transverse projector in the expression of Ji nt . The reason is 
that the definition (|4.6|) restricts anyway the integration to divergenceless fields. 

Finally, we introduce also the bare generating functional which is Q4.31Q when I = I , 

Z (l,m) = -i- [du(x)D P (4.35) 

Zjq J 

X 



exp |z/ Q + J dx[u(x) ■ I(x) + p(x) • m(x)]| 



This bare generating functional is used in the formulation of the perturbative expansion. 
We will first use it to evaluate the bare propagators. 



C. The evaluation of the bare propagators 

It is convenient to change the variables in Jo from p(x) and u(x) to p(k, u) and u(k, u), 
which are defined according to Eq. fl2.19| ). The Jacobian of the transformation (which is 



formally a divergent constant) cancels with the partition sum: 

f dq 



(27T) 4 

1 f dq 



+ g / jAiPa(q)D aP (q) P p(-q) (4.36) 
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where we remind the reader that q = (k, uj). 

The bare propagators are computed from Z (\, m) which in this presentation is written 

as 

Z (l,m) = i- [Du(q)Dp(q) (4.37) 



Z 

x exp {il + J j^jHl) ■ !(-<?) + p(?) • m(-g)]} . 



This functional integral can be computed explicitly. Since we have no mixture of different 
g's in the exponent, the exponential can be computed (on the grid) as the product of 
exponentials. The functional integration is represented as a product of integrals, each one 
for one discrete q and — q: 



Du(q) = J] / du{q)du{-q) 



g>0' 

Computing the resulting Gaussian integrals, and using the expression of the bare Green's 
function (|2.16| ), we end up with 

Z (l,m) (4.38) 

+ G%(q)la(-q)m p (-q)]} . 

We can check that this zero th order generating functional gives the same results as the 
zero th order quantities defined in the direct perturbation theory. In computing the statistical 
quantities from (|4.38|) we use the general property of functional differentiation 

5l a (q) 



8l p {q>) 



(27r) 4 5(g - q')5 a p . (4.39) 



and the symmetry of the bare Green's function G° a Jq) = —G Q *J—q) [cf. Eq. (|2.16| )1. For 



example, defining the velocity correlation function in q- representation, F a p(q), according to 

(u a (q)u p (q')) = (2n) A S(q + q')F a p(q) , (4.40) 
we compute its bare value (u a (q)up(q'))° as 

In other words, we find 
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F%{q) = G ai (q)D jS (q)G* Sf) ( q ) . 



(4.42) 



Similarly, defining the Green's function in g-representation as 



(27r) 4 5(g + q')G aP {q) 

= . s(u a (g)) 

" * Shp(q') 



{u a (q)pp(q')) ■ 



(4.43) 



Since 



(u a (q) P0 (q'))° 



5 2 Z (l,m) 



(4.44) 



8m a (-q)8lp(-q') 



we get trivially that 



(u a (q)Pf3(q')) 



{2ir)% q + q')G%{q). 



We see that this theory generates the same zeroth order propagators as the direct per- 
turbation expansion of Section II. Note however that there exists an apparent additional 
propagator here which is absent in II, which is {p a (q)P/3(q'))- However, the zeroth order 
{Pa(q)Pp(q'))° is zero. This follows from the fact that Z (l, m) does not have a term quadratic 
in m in the exponent. We shall show later that it is zero to all orders. 

Next we show that the statistics inherited from Z (\, m) are Gaussian statistics. The 
exponent contained in Z (l, m) can be expanded as 



Clearly, if we take an odd number of derivatives with respect to l a , and then send 1 and m 
to zero, the result vanishes. On the other hand, if we take 2n derivatives, and then send 1 
and m to zero, only the contributions coming from the n-th power of the integral survive. 
The answer will be proportional to [F°] n with n delta functions. The number of terms will 
be 2n\/[n\2 n ] which is precisely the (2n — 1)!! number of terms expected in the Gaussian 
statistics (cf. Sect. II C). Note that we have Gaussian statistics also for the field p. 




+L(-q)G%(q)mp(-q) 



(4.45) 
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Instead of looking at the correlation functions, one can study the cumulants. These of 
course need to vanish beyond the lowest order. An immediate way to obtain the cumulants 
is to take functional derivatives of In Z (l, m) instead of Z (\, m) itself. Obviously, 

InZ (l,m) = / ^ (4.46) 



l a (-q)F° p ( q )lp(-q) + l a (- q )Gl p (q)me(-q)] . 



Only the second derivative with respect to 1 survives here, giving us only the bare 2-point 
propagators F®p and G ^ as the non-vanishing cumulants, as is required for Gaussian statis- 
tics. 

(a) 

G = \TJV + 1/2 ^ v/W- ^wS' ^ 4- 

X ] 1/2 I 1/2 X ^ T 



4+< x- Yn Y Y Y ~*> 



X] x 2 x 3 x 4 



(b) 

lAAT • 

x x' 
- B 



AG 4 = 




X Xj X2 X 1 

E 



FIG. 12. The graphic representation of Eq.(4.50). Panel (a): The bare, the second order and 
the fourth order terms. Panel (b): All the diagrams up to order F 2 and (c) representative diagrams 
of order T 4 that did not appear in the Wyld diagrammatic series in Fig. 5. All these new diagrams 
vanish. 



AG 2 = 




x x' 
A 
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D. Diagrammatic expansion 



We return to the generating functional Z(\, m) of Eq. (|4.31 ) and write now the propagator 
F a p and G a p as 



F a p(x,x') 
G a p(x,x') 



(u a {x)up{x')) 
(u a (x)up(x') exp{i/ int })° , 
(u a (x)pp(x')} 
(u a (x)p p (x') exp{i/ int })° , 



with il int as defined as [cf. (|2.18| ) and (f4.34|) 1 

Hint = \j dx P^){ r ZU}c 



Expanding the exponential in ( }4.47|) and ( [4.48|) we find 



F = x rjw — + 1/2 ^ \jw- '^J^, C ^J^, ,JVU ti 

x l x 2 



(a) 

> 



+ 4T < <aat- 
x 



jl/2 jl/2 1 1/2 |l/2 x ^ 



x l X 2 x 3 X 4 



(b) 

xi x- 



AF 2 = _ X| r "1 X2 XlC W 2 



X X 

A 



+ 


AAAAAA 




x x' 




" B 



(c) 



AF 4 = 






X 1 ^AA/ V - — -i x 4 






+ 






AAaaan 




x x' 




D 


+ 



(4.47) 
(4.48) 



(4.49) 
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FIG. 13. The graphic representation of Eq.(4.51). Panel (a): The bare, the second order and 
the fourth order terms. Panel (b): All the diagrams up to order T 2 and representative diagrams of 
order T 4 that did not appear in the Wyld diagrammatic series in Fig. 6. All these new diagrams 
vanish. 

G a/3 (X,X') = —( u a{x)pp{x') 

x{\jdxp a ( X )[rZ{:]}J)\ (4.50) 
F a/3 (x, x) = — } {u a {x)pp{x) 

' dx Pa {x)[Vl^J)\ (4.51) 



Such an expansion contains products of the fields u and p averaged over the Gaussian en- 
semble. As explained above, only products of {u a (x)pp(x')) Q and (u a (x)up(x'))° can appear 
in such a series. (Remember that (p a {x)Pi3{x'))° = 0). As done in chapter II, the procedure 
of computing the averages can be represented diagrammatically by pairing all the allowed 
pairs of fields and gluing. The needed graphic notation is the same as the one shown in 
Fig. [I]. The field u is represented by a short wavy line, the field p by a short straight line, 
and now it is obvious why the double correlator is represented by a twice long wavy line. 
Also, the notation for the Green's function becomes natural now, being a correlator of u 
and p. In addition, the notation for the vertex as a dot connecting one straight and two 
wavy tails becomes apparent in Eqs. ( f4.50| ) and fl4.51| ). 

The diagrams that are obtained from Eqs. ( |4.50| ) and ( |4.51| ) are shown in Figs. |T^ a and 



13 a respectively. Upon pairing all the possible pairs and gluing, we get exactly the diagrams 



appearing in Figs. |a and |]a respectively, and in addition we get the diagrams appearing 
in panels b and c of Figs. [12] and [13]. These additional diagrams, that are denoted here 
as AG 2 , AG 4 , AF 2 and AF 4 respectively, do not appear in the diagrammatic expansion 
leading to Figs. || and || However, it is easy to see that they are all zero. The diagram 
A in |I~2] b contains one term which is a diagram for (u) and another one which represents 
a contribution to (p). The latter contains a frequency integral over G°, which is —i/2. 
However the vertex T in front of this contains a factor of k which is odd in the k-integration. 



Therefore the diagram vanishes. The diagrams B, C, D and E in |12|b,c contain closed loops 
of Green's functions. Such loops vanish because of causality: in t-representation only t = 
is allowed for the Green's function in a closed loop, and thus the integral over time has a 
contribution only from a single point, and it vanishes. 
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Next we argue that this phenomenon is general: all the diagrams that do not appear 
in Wyld's expansion are zero. These diagrams are either unlinked (like all the diagrams in 



Fig. 12 except for E) or linked to the body via wavy lines (like E) . All these diagrams must 
have a closed loop of Green's functions. Consider a diagram from the n-th order, in which 
m vertices, m < n, form a fragment that is unlinked to the body. Glue the straight line of 
the first vertex to the wavy line of another. Then either we close the straight line of the 
second vertex on the wavy line, or we connect it to a wavy line of a third vertex. The first 
possibility leads to zero due to the argument of diagram A in |12|b. In the second possibility 
we again need to repeat the same option, or to closed with the first vertex. The production 
of a close loop of Green's functions is inevitable, and these vanish due to causality. Diagrams 
that are connected only via wavy lines also must have a closed loop of Green's functions for 
the same reasons. 

We thus conclude that the path integral formulation leads to exactly the same diagrams 
as the Wyld direct perturbative expansion. This fact is not widely recognized. One reason 
for this is that in Ref. |9| the three types of vertices which appear in this theory were not 



explictly discussed. They were first dealt with properly in [Kj in which it was stated that 
the results agree with || only to order T 3 . We have shown above that the two techniques 
generate exactly the same (non-zero) diagrams on the naive level, and these can be resummed 
in exactly the same way. For example, Fig. |6| of Ref [21] was claimed to differ from the Wyld 
technique, but in fact is identical to the present Fig. |lO| b> which was obtained entirely within 
the Wyld expansion. 



5. SUGGESTIONS FOR FURTHER READING 

The material presented in these notes should be sufficient to prepare the student for 
reading the research papers on the subject of renormalized perturbation theory of turbulence. 
For the convenience of the interested student we present in this section a brief guided tour 
through the latest publications of our own group. The next paper to read is [16|] in which the 



Belinicher-L'vov renormalization scheme is used to prove that the theory for the correlation 
functions and the Green's functions is finite. In other words, it is shown (order by order) 
that for Re— ■> oo, after the appropriate resummations, all the diagrams converge both in 
the IR and the UV limits. This fact tells us that there is no length scale available to 
form dimensionless corrections to the K41 scaling. In |L9[ the issue of anomalous scaling is 



discussed. It is shown that there exist anomalous exponents in turbulence, but they appear 
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in gradient fields and their correlations. Thus for example the dissipation field is anomalous, 
and its correlation function is strongly renormalized compared with the Gaussian limit. The 



issue of the effect of the anomalous fields on the structure functions was discussed in ||20|| . It is 
argued there that the anomalous fields are subcritical, and that they influence the structure 
functions only on the level of corrections to scaling. The latter corrections are large, and 
they disappear only at very large values of Re. We attributed the observed deviations in 
experiments to these "subciritcal" corrections. Two more papers that may be considered 
are [22| and [23] in which further relations to experiments are discussed. In we discuss 
which length scale should be used as the renormalization scale in turbulence, i.e the inner 



scale rj or the outer scale L. In [gz] one can find an experimental test of some aspects of the 
new theory, especially some unusual correlation functions between velocity differences and 
the dissipation field. 
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6. APPENDIX: THE JACOBIAN 



In this appendix we explain the known fact fl7| , |18|l that the Jacobian that appears in 



Eq.(4.10) is unity. From the definition of the delta function it is obvious that 

J Dn(x)6[u(x) -Af-^iix)}] = 1 . (6.1) 

This means that also the functional integral over the RHS of Eq.(4.10) is unity. What we 
want to prove now is that also the functional integral on the RHS of (4.10) without </{u} is 
unity: 

J Du{x)5[Ar{u{x)} - f(x)] = 1 . (6.2) 

This will be the first indication that the Jacobian J{u} = 1. To accomplish this we will 
only use the fact that the Navier-Stokes equations are first order in the time derivative: 

N{u(x)} = + NL{u(x)} . (6.3) 
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The time derivative on the time grid may be represented in the retarded convention 

5^ = [u B (r)-u B _ 1 (r)]/T. (6.4) 
Using this representation we rewrite ( |6.2|) on the time grid as 



/ HDu n (r)5[u n (r) + NLK^r)}] . (6.5) 

We have collected everything except u n itself into an operator that was denoted as NL. It 
is obvious now that every Z)u n (r) integration yields unity due to the delta function, and 
finally the result is unity. Next one should prove that the functional integrals ( |6.1|) and 
( |6.2|) but with any arbitrary weight functional of u are also the same. Repeating the same 
considerations on the time grid furnishes this demonstration with essentially the same ease. 
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